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CUBIC DOMAIN 



HUIYUAN LI, JIACHANG SUN, AND YUAN XU 



Abstract. Several cubature formulas on the cubic domains are derived using 
the discrete Fourier analysis associated with lattice tiling, as developed in ) 101 . 
The main results consist of a new derivation of the Gaussian type cubature for 
the product Chebyshcv weight functions and associated interpolation polyno- 
mials on [— 1,1] 2 , as well as new results on [—1,1] 3 . In particular, compact 
formulas for the fundamental interpolation polynomials are derived, based on 
n 3 /4 + C(n 2 ) nodes of a cubature formula on [—1, l] 3 . 



1. Introduction 

For a given weight function W supported on a set O G M. d , a cubature formula of 
degree 2n — 1 is a finite sum, L n f, that provides an approximation to the integral 
and preserves polynomials of degree up to 2n — 1; that is, 

N 

/ f(x)W(x)dx = V X k f(x k ) =: L n f for all / e 11^, 
Jn fc=i 

where Ilf j denotes the space of polynomials of total degree at most n in d variables. 
The points x k S M. d are called nodes and the numbers € R\{0} are called weights 
of the cubature. 

Our primary interests are Gaussian type cubature, which has minimal or nearer 
minimal number of nodes. For d = 1, it is well known that Gaussian quadrature of 
degree 2n — 1 needs merely N = n nodes and these nodes are precisely the zeros 
of the orthogonal polynomial of degree n with respect to W . The situation for 
d > 1, however, is much more complicated and not well understood in general. 
As in the case of d = 1, it is known that a cubature of degree 2n — 1 needs at 
least N > dimll^-L number of nodes, but few formulas are known to attain this 
lower bound (see, for example, [TUTU]). In fact, for the centrally symmetric weight 
function (symmetric with respect to the origin), it is known that the number of 
nodes, N, of a cubature of degree 2n — 1 in two dimension satisfies the lower bound 

n I 

(1.1) JV>dimn2_ 1+ [-J , 

known as Moller's lower bound [TJJ . It is also known that the nodes of a cubature 
that attains the lower bound (jl.ip . if it exists, are necessarily the common zeros 
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of n + 1 — |_§J orthogonal polynomials of degree n with respect to W . Similar 
statements on the nodes hold for cubature formulas that have number of nodes 
slightly above Moller's lower bound, which we shall call cubature of Gaussian type. 
These definitions also hold in d-dimension, where the lower bound for the number 
of nodes for the centrally symmetric weight function is given in [12] . 

There are, however, only a few examples of such formulas that are explicitly 
constructed and fewer still can be useful for practical computation. The best known 
example is f2 = [—1, l] d with the weight function 

d * d 

(1.2) Woix ) := Y[^=^ or Wi{x) := [] y 1 

and only when d = 2. In this case, several families of Gaussian type cubature are 
explicitly known, they were constructed ([Ml [17]) by studying the common zeros of 
corresponding orthogonal polynomials, which are product Chebyshev polynomials 
of the first kind and the second kind, respectively. Furthermore, interpolation 
polynomials bases on the nodes of these cubature formulas turn out to possess 
several desirable features ([Ml, and also [5]). On the other hand, studying common 
zeros of orthogonal polynomials of several variables is in general notoriously difficult. 
In the case of (|1.2|) , the product Chebyshev polynomials have the simplest structure 
among all orthogonal polynomials, which permits us to study their common zeros 
and construct cubature formulas in the case d = 2, but not yet for the case d = 3 
or higher. 

The purpose of the present paper is to provide a completely different method 
for constructing cubature formulas with respect to Wo and W\ . It uses the discrete 
Fourier analysis associated with lattice tiling, developed recently in [10]. This 
method has been used in [10] to establish cubature for trigonometric functions on 
the regular hexagon and triangle in R 2 , a topic that has been studied in [T51 ITS] , 
and on the rhombic dodecahedron and tetrahedron of R 3 in [9] . The cubature on 
the hexagon can be transformed, by symmetry, to a cubature on the equilateral 
triangle that generates the hexagon by reflection, which can in turn be further 
transformed, by a nontrivial change of variables, to Gaussian cubature formula 
for algebraic polynomials on the domain bounded by Steiner's hypercycloid. The 
theory developed in [TU] uses two lattices, one determines the domain of integral and 
the points that defined the discrete inner product, the other determines the space of 
exponentials or trigonometric functions that are integrated exactly by the cubature. 
In [TO] [9] the two lattices are taken to be the same. In this paper we shall choose 
one as 1 d itself, so that the integral domain is fixed as the cube, while we choose 
the other one differently. In d — 2, we choose the second lattice so that its spectral 
set is a rhombus, which allows us to establish cubature formulas for trigonometric 
functions that are equivalent to Gaussian type cubature formulas for Wq and W\. 
In the case of d = 3, we choose the rhombic dodecahedron as a tiling set and obtain 
a cubature of degree 2n — 1 that uses n 3 /4 + 0(n 2 ) nodes, worse than the expected 
lower bound of n 3 /6 + 0(n 2 ) but far better than the product Gaussian cubature of 
n 3 nodes. This cubature with n 3 /4 + C(n 2 ) nodes has appeared recently and tested 
numerically in [7]. We will further study the Lagrange interpolation based on its 
nodes, for which the first task is to identify the subspace that the interpolation 
polynomials belongs. We will not only identify the interpolation space, but also 
give the compact formulas for the fundamental interpolation polynomials. 
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One immediate question arising from this study is if there exist cubature formulas 
of degree In — 1 with rt 3 /6 + 0(n 2 ) nodes on the cube. Although examples of 
cubature formulas of degree 2n — 1 with N = dirull*^ = n d /d\ + 0(n d ~ 1 ) nodes 
are known to exist for special non-centrally symmetric regions ([1]), we are not aware 
of any examples for symmetric domains that use N = n d jd\ + 0(n d ~ 1 ) nodes. From 
our approach of tiling and discrete Fourier analysis, it appears that the rhombic 
dodecahedron gives the smallest number of nodes among all other fundamental 
domains that tile R 3 by translation. Giving the fact that this approach yields the 
cubature formulas with optimal order for the number of nodes, it is tempting to 
make the conjecture that a cubature formula of degree 2n — 1 on [—1, l] 3 needs at 
least n 3 /4 + 0(n 2 ) nodes. 

The paper is organized as follows. In the following section we recall the result 
on discrete Fourier analysis and lattice tiling in [lOj . Cubature and interpolation 
for d — 2 are developed in Section 3 and those for d = 3 are discussed in Section 4, 
both the latter two sections are divided into several subsections. 

2. Discrete Fourier Analysis with lattice Tiling 

We recall basic results in [10] on the discrete Fourier analysis associated with a 
lattice. A lattice of R d is a discrete subgroup that can be written as AZ d = {Ak : 
k G Z d }, where A is a d x d invertible matrix, called the generator of the lattice. A 
bounded set n A C R d is said to tile R d with the lattice AZ d if 

xn A (x + Ak) = 1 for almost all x G K d , 

where \ E denotes the characteristic function of the set E. The simplest lattice is 
Z d itself, for which the set that tiles M. d is 

fl:= [-±,±) d . 

We reserve the notation as above throughout the rest of this paper. The set f2 
is chosen as half open so that its translations by Z d tile M. d without overlapping. It 
is well known that the exponential functions 

e fc (a:) :=e 2nik - x , k G Z d , xeR d , 

form an orthonormal basis for L 2 (fl). These functions are periodic with respect to 
Z d ; that is, they satisfy 

f(x + k) = f(x) for all k G Z d . 
Let B be a d x d matrix such that all entries of B are integers. Denote 
(2.1) K B = {k G Z d : B~ tr k G ft} and A^ = {k G Z d : k G fl B } ■ 

It is known that |As| = |A^| = |det-B|, where \E\ denotes the cardinality of the 
set E. We need the following theorem [10l Theorem 2.5]. 

Theorem 2.1. Let B be a d x d matrix with integer entries. Define the discrete 
inner product 

{f,g)B ■■= E f{B-«jWB^T) 
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for f, g £ C(Q), the space of continuous functions on fl. Then 

(2.2) </, g) B = (/, <?) := / f(x)gjx)dx, 

Jn 

for all f, g in the finite dimensional subspace 

T B := span { e 2mkx : k£ A^} . 

The dimension of Tb is |A^| = |det£>|. 

This result is a special case of a general result in [TU] , in which f2 is replaced by 
Ha for an invertible matrix A, and the set A^ is replaced by Km with N = B tr A 
and N is assumed to have integer entries. Since we are interested only at the cube 
[— i, i] d in this paper, we have chosen A as the identity matrix. 

We can also use the discrete Fourier analysis to study interpolation based on the 
points in As. We say two points x,y £ M d congruent with respect to the lattice 
BJj d , if x — y £ BZ d , and we write x = y mod B. We then have the following 
result: 

Theorem 2.2. For a generic function f defined in C(Q), the unique interpolation 
function Tb f in Tb that satisfies 

I B f(B- tr j) = f(B^j), Vj £ A B 

is given by 

(2.3) T B f(x) = ]T (f,e k )e k (x) = £ /(£r tr fc)*n B (x - B^k), 

fceA]g fe eA B 

where 

The proof of this result is based on the second one of the following two relations 
that are of independent interests: 

(25) 1 c 2*ik u B- v i _ [l> iffc = modB, 

l ' J \det{B)\ ^ \Q, otherwise, 

and 

/ 2 ^ 1 e -27rife»s-*'j = f 1, if J = mod £ tr , 

|det(S)| ^ |0, otherwise. 

fceA T „ *> 



For proofs and further results we refer to [101 [9] . Throughout this paper we will 
write, for k £ Z d , 2k = (2fci, . . . , 2k d ) and 2k + 1 = (2ki + l,...,2k d + 1). 

3. CUBATURE AND INTERPOLATION ON THE SQUARE 

In this section we consider the case d = 2. In the first subsection, the general 
results in the previous section is specialized to a special case and cubature formulas 
are derived for a class of trigonometric functions. These results are converted to 
results for algebraic polynomials in the second subsection. Results on polynomial 
interpolation are derived in the third subsection. 
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3.1. Discrete Fourier analysis and cubature formulas on the plane. We 

choose the matrix B as 



B = n 



1 1 

-1 1 



and 



B- 1 = 



1 

2n 



1 -1 
1 1 



Since B is a rotation, by 45 degree, of a constant multiple of the diagonal matrix, 
it is easy to see that the domain fig is defined by 

Sis = {x £ I 2 : — n < x\ + x 2 < n, —n < X2 — x\ < n}, 

which is depicted in Figure 1 below. 



Figure 1. Rhombus fin 

From the expression of B~ tr , it follows readily that Ab = A B =: A„, where 

A„ = {j e Z 2 : -n < ji + j 2 < n, -n < j 2 - ji < n}. 

The cardinality of A n is |A n | = 2n 2 . We further denote the space Tb by T n , which 
is given by 

T n :=span{e 2 " fe ;r : k € A„} . 
Theorem 3.1. Define the set 

X n := {2k : -f < k u k 2 < f } U {2k + 1 : -2±i < k x ,k 2 < ^} . 
Then for all f,g eT„, 

Proof. Changing variables from j to k = 2nB~ tr j, or k± = ji + j 2 and k 2 = j 2 — j\, 
then, as j\ and j2 need to be integers and j\ — kl ~ k2 , j 2 = kl + k2 ; we see that 

(3.1) j e A n k = 2nB- tr j e X n . 

Hence, as det(_B) = 2n 2 , we conclude that (/, g) n — (/, g) B and this theorem follows 
as a special case of Theorem 12.11 □ 

The set A„ lacks symmetry as the inequalities in its definition are half open and 
half closed. We denote its symmetric counterpart by A* , which is defined by 

A* := {j £ I? : -n < ji + j 2 < n, -n < ji - j 2 < n}. 
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We also denote the counterpart of T n by T* , which is defined by 

T„*:=span{e 27r4fca; : fc <E A*} . 

Along the same line, we also define the counterpart of X n as 

X* := {2k : -f < fci, k 2 < § } U {2k + 1 : -2±i < fci, k 2 < ^}. 

It is easy to see that \X n \ = |A„| = 2n 2 , whereas \X*\ = 2n 2 + 2n + 1. We further 
partition the set X* into three parts, 

x n = x° U X!„ u x v n , 

where X° = X* n (— n, ra) 2 is the set of interior points of X*, X% consists of those 
points in X* that are on the edges of [— n, n] 2 but not on the 4 vertices or corners, 
while X% consists of those points of X* at the vertices of [—n, n] 2 . 



Theorem 3.2. Define the inner product 

(3-2) if, 9)1 4l 4"V(!rMl), when 
kex* 

Then for all /, g €iT n , 



c 



1, kex° 

2 , K t ^„ . 
41 ^ ^ -^n 



f{x)g{x)dx = (f,g) n = (f,g)* n . 



Proof. Evidently we only need to show that (f,g) n = (f,g)n- Since cj^ = 1 for 
k e X°, the partial sums over interior points of the two sums agree. The set 
X^ of boundary points can be divided into two parts, X^ = X^ 1 U X^ 2 , where 
X^ 1 consists of points in X n that are on the edges of [— n,n) 2 , but not equal to 
(— n, — n), and X^ 2 is the complementary of X^ 1 in X^. Evidently, if x € X^ 1 , 
then either x + (2ra, 0) or x + (0,2n) belongs to X^ 2 . Hence, if / is a periodic 
function, f(x + k) = f(x) for k e Z 2 , then 

E <£°/(&) = 5 E /(&) - E /(&)■ 

kex- keX e,i 

Furthermore, for (— n, — n) € X n , X* contains all four vertices (±n, ±n). Since 
a periodic function takes the same value on all four points, J2kex^ c fe"' > /(^r) = 
/(— 5, —5)- Consequently, we have proved that (/, g) n — (/, g)* n if /, g are periodic 
functions. □ 

As a consequence of the above two theorems, we deduce the following two cuba- 
ture formulas: 

Theorem 3.3. For n>2, the cubature formulas 
(3.3) 

and / 11]2 ^ = iE/^) 

"'L-s-sJ fcex* ■ / 1-2>3J keXn 

are exact for f £ T^n-i- 
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Proof. It suffices to proof that both cubature formulas in (|3.3p are exact for every 
ej with j E A|„_i. For this purpose, we first claim that for any j 6 Z 2 , there exist 
i/£ A„ and Z g Z 2 such that j — v + Bl. Indeed, the translations of by -BZ 2 
tile IR 2 , thus we have j = x + Bl for certain x E and I E Z 2 . Since all entries of 
the matrix B are integers, we further deduce that v := x = j — Bl E I? <1£Ib = A n . 

Next assume j E A2 n _ 1 . Clearly the integral of e^- over fl is Sjfl. On the other 
hand, let us suppose j — v + Bl with v € A„ and i £ Z 2 . Then it is easy to see that 
e j(^l) = ev (w~) ^ or eac ^ ^ *= ^n- Consequently, we obtain from Theorem 13.21 that 

fcex* fcex* fcex„ 

= E e j(^") = / e„(x)da; = ^ l0 . 

Since z/ = implies j = Bl E Z 2 which gives j = I = 0, we further obtain that 
$v,o = Sj.o- This completes the proof of (|3.3p . □ 

We note that the second cubature in (|3.3[) is a so-called Chebyshev cubature; 
that is, all its weights are equal. 

3.2. Cubature for algebraic polynomials. The set A* is symmetric with re- 
spect to the mappings (a; 1,0:2) l— * {—x\,X2) and (351,0:2) (x±, — #2)- It follows 
that both the spaces 

7^ ven : = span{cos27rjia;i cos27rj2X2 : < ji + j'2 < n}, 

T° dd : = span{sin27rjixi sin27rj2^2 : 1 < ji + 32 < n} 

are subspaces of 7^*. Recall that Chebyshev polynomials of the first kind, T n (t), 
and the second kind, U n (t), are defined, respectively, by 

TJt) = cosnO and UJt) = &m ( n + jOg t = cos6 ,_ 

sin ft 



They are orthogonal with respect to w (t) = 1 / >/l — t 2 and iui(t) = \J\—t 2 over 
[—1,1], respectively. Both are algebraic polynomials of degree n in t. Recall the 
definition of Wo and W\ in (|1.2|) . Under the changing of variables 

(3.4) t 1 = cos27rxi, t 2 — cos 27ra 2 , (xi, x 2 ) € [— \, ^] 2 , 

the subspace 7^ vcn becomes the space II 2 of polynomials of degree n in the variables 
(ti.ta), 

n 2 = span^^tiJTfc-j-Cta) : < j < k < n} 

and the orthogonality of et over f2 implies that TL- (i) := Tj(ti)Tk-j(t2) are orthog- 
onal polynomials of two variables, 

fl, fe = fc'=i = i' = 0, 
1 /" , , i = (fc'.i') and (k-j)j = 0, 

^ J[-isy \ jj k = k > j = j > Q, 



We note also that the subspace T° dd becomes the space {yl — t\\J\ — t\-pit) : p E 
H-n-i} m the variables t = (ii, £2)3 and the orthogonality of efe also implies that 
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UHt) := Uj(ti)Uk-j{t2) are orthogonal polynomials of two variables, 



' 1 tf{t)Uf,(t)W 1 (t)dt=- A 5 jtj ,5 k ,k> 
[-M] 2 4 



The symmetry allows us to translate the results in the previous subsection to 
algebraic polynomials. Since cos 2-njiX\ cos 2irj 2 x 2 are even in both variables, we 
only need to consider their values over X* n {x : X\ > 0, x 2 > 0}. Hence, we define 

(3.5) 3„ := {(2fci, 2fea) : < k u k 2 < f }U{(2fci + 1, 2fc 2 + 1) : < k u k 2 < 



and, under the change of variables (|3.4p . 

(3.6) T„ := {(z kl ,z k2 ) : (h,k 2 ) € S„}, where z fe = cos^. 

Furthermore, we denote by T° := r„ n (—1, l) 2 the subset of interior points of r„, 
by the set of points in T n that are on the boundary of [— 1, l] 2 but not on the 
four corners, and by the set of points in r„ that are on the corners of [— 1, l] 2 . 
The sets 3°, 5^ and 5^ are defined accordingly. A simple counting shows that 



.§'•>• 



(3.7) |5„| = ([|J + l) 2 + ([^j+l)' = lfctil 
Theorem 3.4. The cubature formula 

(3.8) \ f f(t)W (t)dt=^ **°/(*kx,*fe,), 
is exact for II| n _ 1 . 

Proof. We note that X* is symmetric in the sense that k £ X* implies that 
(— ki,k 2 ) € X* and (ki,—k 2 ) € X*. Let g(x) = /(cos 2-kxx, cos 2^X2). Then 3 
is even in each of its variables and g(-Mr) = f(z kl ,Zk 2 ). Notice that / € n| n _ 1 
implies g € T 2n -i- Applying the first cubature formula (|3.3p to <?(x), we see that 
follows from the following identity, 




To prove this identity, let ka denote the set of distinct elements in {(±ki,±k 2 )}; 
then 5(5^) takes the same value on all points in ka. If k S X*, fci ^ and fc 2 7^ 0, 



2g{±) if k e X«, and £. £fe(T c< n) ff (£) = if k E X*. If h = and fc 2 ^ or 

fc 2 = and fci 7^ 0, then ka contains 2 points; J2jek<j c i™' ^ d(4rd ~ ^Q^lrd ^ ^ e "^n 
and Eieto c k9{4^) = 9{^) if S X*. Finally, if fc = (0,0) then fc<r contains 1 
point and 5(0, 0) has coefficient 1. Putting these together proves the identity. □ 



By (|3.7p , the number of nodes of the cubature formula (|3 . 8|) is just one more than 
the lower bound p.l[) . We can also write (|3.8p into a form that is more explicit. 
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Indeed, if n — 2m, then (13. 8p can be written as 



(3.9) -1 f f(t)W (t)dt 
n J [-1,1]' 



mm m _i m _i 

2 v-^" v-^" , x 2 



where J^" means that the first and the last terms in the summation are halved. If 
n = 2m + 1 , then (13.81) can be written as 



(3.10) \ I f(t)W Q (t)dt 



mm mm 

2 W W „, % 2 



W * — ' * — ' n" 

where ^ means that the first term in the sum is divided by 2. The formula (|3.10p 
appeared in [17j . where it was constructed by considering the common zeros of 
orthogonal polynomials of two variables. 

From the cubature formula (|3.3|) . we can also derive cubature formulas for the 
Chebyshev weight W\ of the second kind. 

Theorem 3.5. The cubature formula 



- 2 \ f(t)w 1 (t)dt=^ sii 



(3.H) — / f(t)Wx(t)dt= — ^ si n^sin^/(z fcl ,z fe2 

fc£H° 

is eiaci /or Il| n _ 5 . 

Proof. We apply the first cubature formula in (|3 . 3|) on the functions 

sin(27r(fci + l)xi) sin(27r(fc2 + 1)2:2) sin 2ttx\ sin 2itx2 

for < fci + &2 2n — 5, where ii = cos27ra;i and £2 = cos27TX2 as in 
Clearly these functions are even in both x\ and x 2 and they are functions in T^n-i- 
Furthermore, they are zero when x\ — or x 2 — 0, or when (x\,X2) are on the 
boundary of X*. Hence, the change of variables (|3 .4|) shows that the first cubature 
in $£5$ becomes (|3~TT|) for U kl (ti)U k2 (t 2 ). □ 

A simple counting shows that |S°| = LfJ 2 + L^J 2 = ( "~ 1) 2 ( "~ 2) + [|J. The 
number of nodes of the cubature formula (|3.1ip is also one more than the lower 
bound (|1.1|) . In this case, this formula appeared already in [13]. 

3.3. Interpolation by polynomials. As shown in [10] , there is a close relation 
between interpolation and discrete Fourier transform. We start with a simple result 
on interpolation by trigonometric functions in T n . 

Proposition 3.6. For n > 1 define 

(3.12) I n f(x) := - * n (x) ^2 E e ^ 

Then Inf(^) = /(£) for all k e X n . 
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Proof. For j G A„ define k = 2nB tr j. From the relation (|3.1I) . j € A n is equivalent 
to k G X n with fc = 2nB~ tr j. As a result, we can write I n f(x) as 

j'eA„ 

and the interpolation means Inf(B~ tr j) = f(B~ tr j) for j G A„. For j G A„, 



$„(S- tr (j - k)) = -L £ e,(5- tr (j - fc)) 



by (EH). □ 

For our main result, we need a lemma on the symmetric set X* and A* . Recall 
that cjf 1 is defined for k G X*. Since the relation (|3.ip clearly extends to 

(3.13) jeAJ ^ fc = 2nB~ tr j G A"*, 

we define cj- = whenever fc and j are so related. Comparing to (I3.12| . we 
then define 

(3.14) I*J{x) := J2 f&li* ~ wh ere $;( x) := _i_ ^ ^(x). 

fcex* ^eA* 

We also introduce the following notation: for fc G A"^, we denote by k 1 the point on 
the opposite edge of X*; that is, k' G X e n and k' = k ± (2n, 0) or fc' = fc ± (0, 2n). 
Furthermore, we denote by j' the index corresponding to k' under (|3 . 1 3(1 . 

Lemma 3.7. The function I*f G T* satisfies 



I* f(JS-) = I 



/(£)> k G XI 



/(£)+/(£)> 



U(li) + /(^ M ) + /(^ M ) + /(^) I *e**. 

Proof. As in the proof of the previous theorem, we can write I*f as 

E /(^- tr i)$;^ - s- tr j) 

by using (j3~T3)l . Let 5* (at) = $*(S- tr j). For all fc, j G A;, 

S k (B- tr j) = ^ 4 n) e,(S- tr (j - fc)). 
i/eA* 

Since e v (B~ tr j) — en{B~ tr j) for any [i = v mod _B, we derive by using a similar 
argument as in Theorem 13.21 that 



Sk (B-«j) = ^ E e u (B^(j-k)). 



2n 2 

^eA„ 



By dHJ), S k (B- tr j) = 5 k -j if fc,j G A n . If j G A* \ A„ then j' G A„, so that 
if k G A„ then Sk{B~ tr j) = 8k. j' ■ The same holds for the case of j G A„ and 
k G A* \ A n . If both k, j G A* \ A„, then S k {B^ tr j) = 5 k >, r . Using the relation 
(|3.13[) , we have shown that $^(2^) = 1 when k = j mod 2nZ 2 and otherwise, 
from which the stated result follows. □ 
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It turns out that the function <&* satisfies a compact formula. Let us define an 
operator V by 

{^f){ x ) = \ [f(xi,x 2 ) + f{-x x ,x 2 ) + f{x\,-x 2 ) + f(-Xi,-X 2 )] ■ 



it follows immediately that 



For e k {x) 

(3.15) (Pek)(x) = cos(27rfcia;i) cos(27rfc 2 ^2) 

Lemma 3.8. For n > 0. 

1 



forall k £ Z 2 . 



(3.16) 

where 
(3.17) 

D n {x) :-- 



$*(x) = 2 [D n (x) + D n -i(x)] — -(cos27rna;i + cos27ma;2), 



4 X! e -( x ) 



1 cos 7r(2n + l)xi cos ixx\ — cos 7r(2n + l)a;2 cos 7nr2 
2 



cos 2irx\ — cos 27ra2 



Proof. Using the values of ci™' ) and the definition of D n , it is easy to see that 



$* n (x)=2[D n (x)+D n _ 1 (x)] 



e v (x). 



Since contains four terms, (±n,0) and (0,±n), the sum over becomes the 
second term in (|3.16p . On the other hand, using the symmetry of A* and l|3.15p . 



i/6A* 0<ji+j 2 <n 



cos 2njiXi cos 2ttj' 2 x 2 , 



where means that the terms in the sum are halved whenever either j\ = or 
j 2 = 0, from which the second equal sign in (|3.17p follows from [18j (4.2.1) and 
(4.2.7)]. □ 

Our main result in this section is interpolation over points in {^- : k £ S„} with 
S ra defined in (|3.5[) . 

Theorem 3.9. For n > define 

C n f(x) = J2 4(x) := A< n >7> [$;(• - A)] (s) 



fces r . 



im'f/i A^ n given in (|3.8p . T/ien £ n / £ T n is even in both variables and it satisfies 

£ «/(^) = /(^r) ies„. 

Proof. As shown in the proof of Proposition 13.71 Rk(x) :— $n(x ~ J~) satisfies 
Rk{4-) = 1 when k = j mod 2nl? and otherwise. Hence, if j £ S° then 

C^fcX^) = 3^(^).= [4" ) ]~ 1 *fcj- If J e then the number of terms 
in the sum of {'PRk)(- 2 ^) depends on whether j\j 2 is zero; if j\j 2 ^ then 

(P*k)(£) = 3 [^C^ + ^fe] = ¥k,j = l^Y^kj, whereas if jij 2 = 

then (VR k )(£) = \Rk{it) = [A* ] -1 ^- For j = (n,0) or (0, n) in S£, we 



have (VR k )(£) 



for j = (n, n) e we have 



Li = I 



_i? fe (^) + ii fe (^) + ii fc (^) + J R fc ((^)] = 4j; fi _ 

nally for j = £ H^, it is evident that (VRk)(0) = 5k,a- Putting these together, we 
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have verified that (-k(o^) = Sk j for all j, k S H* , which verifies the interpolation of 
C n f. '" ' ' □ 

As in the case of cubature, we can translate the above theorem to interpolation by 
algebraic polynomials by applying the change of variables (|3.4p . Recall r„ defined 
in (EH). 



Theorem 3.10. For n>0, let 

£»/(*)= E /(**)**(*). t k{t)=h{x) with U = cos2-kx u i = 1,2. 

T/ien £«/ € 11^ and it satisfies C n f(zk) — f(zk) for all z k E T n . Furthermore, 
under the change of variables (|3 .4[) . the fundamental polynomial l* k {t) satisfies 

m - \V [D n (- - A) + £»„_ l( . _ £)] {x) - i [(-l)*!^) + (-l) fc2 T fc2 (t 2 )] . 

Proof. That £ n / interpolates at € T„ is an immediate consequence of the change 
of variables, which also shows that C n f € n^. Moreover, cos27rn(a;i — |^-) = 
(— l) fel cos27rnxi = (— l) fcl T„(xi), which verifies the formula of t k (t). □ 

The polynomial C n f belongs, in fact, to a subspace n* C IT^ of dimension 
|5 n | = dimll^j + LfJ + 1) and it is the unique interpolation polynomial in II* . In 
the case of n is odd, this interpolation polynomial was defined and studied in [15] , 
where a slightly different scheme with one point less was studied in the case of even 
n. Recently the interpolation polynomials in [19] have been tested and studied 
numerically in [3l |4j; the results show that these polynomials can be evaluated 
efficiently and provide valuable tools for numerical computation. 

4. Cubature and Interpolation on the cube 

For d — 2, the choice of our spectral set and lattice in the previous section 
ensures that we end up with a space close to the polynomial subspace IT^; indeed, 
monomials in 11^ are indexed by < j\ + ji < n, a quarter of A* . For d = 3, 
the same consideration indicates that we should choose the spectral set as the 
octahedron {x : —n < x\ ± X2 ± £3 < n}. The octahedron, however, does not tile 
M. 3 by lattice translation (see, for example, [6] p. 452]). As an alternative, we choose 
the spectral set as rhombic dodecahedron, which tiles K 3 by lattice translation with 
face centered cubic (fee) lattice. In [}J], a discrete Fourier analysis on the rhombic 
dodecahedron is developed and used to study cubature and interpolation on the 
rhombic dodecahedron, which also leads to results on tetrahedron. In contrast, our 
results will be established on the cube [— |, ^] 3 , but our set £Ib is chosen to be a 
rhombic dodecahedron. 

4.1. Discrete Fourier analysis and cubature formula on the cube. We 

choose our matrix B as the generator matrix of fee lattice, 

'0 1 l\ . f-1 1 1 

B = ;, I 1 1 and B- 1 = — 1 -1 1 
T 1 0/ \ 1 1 -1, 

The spectral set of the fee lattice is the rhombic dodecahedron (see Figure 2). Thus, 

B = {x e K 3 : —n < x v ± x^ < n, 1 < v < fi < 3}. 
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(0,0,1) 




(0. 0, -1) 



Figure 2. Rhombic dodecahedron 

The strict inequality in the definition of fl B reflects our requirement that the 
tiling of the spectral set has no overlapping. From the expression of B~ tr , it follows 
that A B — '■ A n is given by 

A„ := {j e I? : -n < -j x + j 2 + j 3 ,ji - h + h, h + h - h < n}. 

It is known that |A„| = det(B) = 2n 3 . Furthermore, =: Aj, is given by 

A^ = Z 3 n fl B = {k G I 3 : -n < k v ± fc M < n, 1 < v < fi < 3}. 

We denote the space T B by T n , which is given by 

T n := span{e 2 " fc :E : k G A*}. 

Then dim 7; = |A+| = det(B) = 2n 3 . 

Theorem 4.1. Define the set 

X n := {2k : -| < k u ka, k 3 < § } U {2k + 1 : -s±i < h, k 2 , k 3 < ^} . 

Then for all /, g G T n , 

Proof. Changing variables from j to k — 2nB~ tr j, or j = B tr k/(2n), then, as 
31,32,33 are integers and j x = ^±^, j 2 = j 3 = k ^ 2 ds2 -, we see that 

(4.1) j G A„ <=► 2nB" tr j G A„ and £ /(£T tr j) = £ /(£), 

jeA„ fcex„ 

from which we conclude that (f,g) n — (f,g) B - Consequently, this theorem is a 
special case of Theorem 12.11 □ 

Just like the case of d = 2, we denote the symmetric counterpart of X n by X* 
which is defined by 

X* := {2k : -f < kx, fc 2 , fc 3 < f } U {2k + 1 : -^±1 < fc 2) k 3 < Sfi}. 

A simple counting shows that \X*\ = n 3 + (n+l) 3 . The set X* is further partitioned 
into four parts, 

x* = x°uxfux°uxi 
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where X° = X* n (— n, n) 2 is the set of interior points, X T { contains the points in 
X* that are on the faces of [— n,n} 3 but not on the edges or vertices, X^ contains 
the points in X* that are on the edges of [— n, n] 3 but not on the corners or vertices, 
while X% denotes the points of X* at the vertices of [— n, n] 3 . 

Theorem 4.2. Define the inner product 



(4-2) </,*>;:= E 4 B) /(&M&), ^ere g> 
kex* 



fl, 


k G X° 


1 

< ?' 


fee A,{ 


4' 


fe G X« 


1 

v 8 ' 


fe G X- 



T/ien /or a/Z /, g G T n , 



i 113 



f{x)g{x)dx = (f,g) n = (f,g)* 



Proof. The proof follows along the same line as the proof of Theorem l3.21 We only 
need to show (f,g) n = (/,<?) if fg is periodic. The interior points of X n and X* 
are the same, so that cj^ = 1 for k € X°. Let e% — (1,0,0), e-i = (0,1,0), and 
£3 = (0, 0, 1). Each point k in X^ has exactly one opposite point k* in X£ under 
translation by ±nei and only one of them is in X n , so that f(xk) = \[f{xk) + f{x* k )] 

if / is periodic, which is why we define c~ k = ^ for k G X£. Evidently, only three 
edges of X* are in X* \ X n . Each point in X* corresponds to exactly four points 
in X^ under integer translations ±nei and only one among the four is in X n , so 
we define = 7 for fc G X%. Finally, all eight corner points can be derived 
from translations nSi points, used repeatedly, and exactly one, (— n, —n, —n), is in 
X* \ X n , so that we define c k n) = f for k G X%. □ 

We also denote the symmetric counterpart of A^ by A^* , 

(4.3) At* := {j G Z 3 : -n < j v ± j„ < n, 1 < v < n < 3} 

and denote the counterpart of T n by T* , which is defined accordingly by 

T„* :=span{e 2 " fc - : k G Aj>* } . 

Theorem 4.3. For n > 2, £/ie cubature formulas 

(4.4) ' 

/ - i E <£°/(&) and / /(*)*: = ^ £ fit) 

are eiaci /or / G T^n-i- 

Proof. As in the proof of Theorem 13.31 for any j G Z 3 , there exist G A^ and 
I G Z 3 such that j = v + Bl. 

Assume now j G A\* n _ 1 . Clearly the integral of e,- over is Sjfi. On the other 
hand, let us suppose j — 1/ + Bl with v G A„ and Z G Z 3 . Then it is easy to see that 
e j(^l) = e,y (^r) f° r eacn ^ € A*. Consequently, we get from Theorem 14.21 that 

E 4 n) e,-(£) = E 4 n) e,(|r) = E e »(£) 

fceA* fcex* feeA„ 



= E e j(^r) = / e »( x )dx = S ufi . 

1, r- v •'a* 



fceA„ 
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Since v = implies j = I = 0, we further obtain that 5 Vt o = Sj^q. This states that 
the cubature formulas (|4.4|) are exact for each e.,- with j € Al* l5 which completes 
the proof. □ 

4.2. Cubature formula for algebraic polynomials. We can also translate the 
cubature in Theorem 14.31 into one for algebraic polynomials. For this we use the 
change of variables 

(4.5) ii = cos2tt:ei, t 2 — cos 2irx 2 , t 3 = cos 27ra;3 , i€[-yf. 

Under (|4.5[) . the functions cos27rfcia;i cos27rfc 2 X2 cos27rfc3a:3 become algebraic poly- 
nomials T kl (ti)Tk 2 {t 2 )T kz (t 3 ), which are even in each of its variables. The subspace 
of T* that consists of functions that are even in each of its variables corresponds 
to the polynomial subspace 

II* := spa,n{T kl (xi)T k2 (x 2 )T k3 (x 3 ) : fa, k 2 ,k 3 > 0, k v + < n, 1 < v < /i < n}. 

Notice that X* is symmetric in the sense that if x G X* then ax e X* for all 
a G { — 1,1} , where (ax)i — OiXi. In order to evaluate functions that are even 
in each of its variables on X* we only need to consider X* n {x : x\,x 2 ,x 3 > 0}. 
Hence, we define, 

(4.6) S n := {2k : < fa, k 2 , k 3 < §} U {2k + 1 : < fa, k 3 < ^} 
and, under the change of variables (|4.5|) . define 

(4.7) r„ := {(z kl ,z k2 ,z k3 ) : k S S„}, z fe = ^. 

Moreover, we denote by T°, r^, and the subsets of T n that contains interior 
points, points on the faces but not on the edges, points on the edges but not on the 
vertices, and points on the vertices, of [— 1, l] 3 , respectively, and we define 5°, S£, 
3® and SJ^ accordingly. A simple counting shows that 

{ (n+l) 3 , 3(ra+l) 
+ 4 , n is even, 
- — n is odd. 

Theorem 4.4. Write z k = (z kl , Zu 2 , Zk 3 )- The cubature formula 



(4.9) 4/ /(^»p4e4"/w. t } 

77 J\-l,l] 3 Ln ke ~ 





fee 


■^7° 
1 — '71 ' 




fee 


a/ 

'—'71' 


2, 


fee 


'—'71' 


.1, 


fee 


'—'71' 



«s exact for Yl^n—i- In particular, it is exact for Il2 n _ 1 . 

Proof. Let g(a;) = /(cos 27rxi, cos 27rx2, cos 2^x3). Then g is even in each of its 
variables and g(^) = f{z k ). Applying the first cubature formula in (|3.3p to g(x), 
we see that (|3.8[) follows from the following identity, 

E 4 n) ff(lr) = E 4 n) iW 

This identity is proved in the same way that the corresponding identity in Theorem 
I3.4l is proved. Let ka denote the set of distinct elements in {kcr : a G {— 1, 1} 3 }; then 
g{^) takes the same value on all points in ka. If k 6 X*, fa 7^ for i = 1,2, 3, then 
ka contains 8 points; if exactly one ki is zero then ka contains 4 points; if exactly 
two fa are zero then ka contains one point; and, finally, if k = (0,0,0) then ka 
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contains one point. In the case of fcj 7^ for i = 1,2, 3, J2j €ktT c k 9{~in) ~ ^(ife) 

if * G E iefcff 4 n) 9(^) = 4 ff(|r) if * e X/, and £ . £fc(T 4" } 9 (^) = 2 ff ( A) if 
G The other cases are treated similarly. Thus, (|4.9[) holds for njj n _i. 
Finally, the definition of IT* shows readily that it contains 

II 3 = sptm{T kl (xi)T k2 (x 2 )T k3 (x 3 ) : h,k 2 , k 3 > 0, < h + k 2 + k 3 <n} 

as a subspace. In particular, H^n-i contains Il2 n _i as a subset. □ 

We note that II| n _ 1 contains Dln-i as a subspace, but it does not contain IT 3 n 
since T n {x\)T n {x 2 ) is in El| n but not in Il2 n _ 1 . Hence, the cubature (|4.9|) is of 
degree 2n — 1. A trivial cubature formula of degree 2n — 1 for Wq can be derived 
by taking the product of Gaussian quadrature of degree 2n — 1 in one variable, 
which has exactly n 3 nodes. In contrast, according to (|4.8[) . the number of nodes of 
our cubature (|3 . 8|) is in the order of n 3 /4 + 0(n 2 ), about a quarter of the product 
formula. As far as we know, this is the best that is available at the present time. 
On the other hand, the lower bound for the number of nodes states that a cubature 
formula of degree 2n—l needs at least n 3 /6 + 0(n 2 ) nodes. It is, however, an open 
question if there exist formulas with number of nodes attaining this theoretic lower 
bound. 

Recall the cubature (|4.9p is derived by choosing the spectral set as a rhombic 
dodecahedron. One natural question is how to choose a spectral set that tiles M 3 
by translation so that the resulted cubature formula is of degree 2n — 1 and has the 
smallest number of nodes possible. Among the regular lattice tiling, the rhombic 
dodecahedron appears to lead to the smallest number of nodes. 

Just as Theorem 13.51 we can also derive a cubature formula of degree 2n — 5 for 
W\ from Theorem 14.31 We omit the proof as it follows exactly as in Theorem 13.51 



Theorem 4.5. The cubature formula 

(4.10) ^ / f(t)Wi(t)dt = 4 E sin2 sin2 *f sin2 ^/(**) 
is exact for II*, n _ 5 ; in particular, it is exact for Il2 n _ 5 . 

4.3. A compact formula for a partial sum. In order to obtain the compact 
formula for the interpolation function, we follow [9] and use homogeneous coordi- 
nates and embed the rhombic dodecahedron into the plane t\ + t 2 + t 3 + t± = 
of R 4 . Throughout the rest of this paper, we adopt the convention of using bold 
letters, such as t, to denote the points in the space 

R% := {t = {t 1 ,t 2 ,t 3 ,t 4 ) S M 4 : h + t 2 + t 3 + t 4 = 0} . 

In other words, the bold letters such as t and k will always mean homogeneous 
coordinates. The transformation between i£l 3 and t G M.jj is defined by 

ii = +x 2 + x 3 ) 

xi=t 2 +t 3 { 

(4.11) {x 2 =t 1+ t 3 — )t 2 = - 2 {x 1 -x 2 +x 3 ) 



x 3 =t 2 +tl 



h = \{x\+x 2 - x 3 ) 
ti = - x 2 - x 3 ). 

In this homogenous coordinates, the spectral set fis becomes 

(4.12) Q B = {t € R% : -1< U - tj < 1, 1 < i < j < 4} . 
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We now use homogeneous coordinates to describe Ajj* denned in (|4.3|) . Let 

Z% := Z 4 n Rjj and 

H := {j e 4 : ji = j 2 = ja = k mod 4}. 
In order to keep the elements as integers, we make the change of variables 

h = 2(-fci + k 2 + k 3 ), 32 = 2(h - k 2 + k 3 ), 
' ' L>l 33 =2{ki + k 2 -k 3 ), ji = 2(-fci - fe 2 - k 3 ) 

for k = (ki,k2,k 3 ) € Ajj*. It then follows that A^* in homogeneous coordinates 
becomes 

G„ := {j E H : ji = j 2 = j 3 = j 4 = mod 2, -4n < j v - < An, 1 < v, \i < 4}. 

We could have changed variables without the factor 2, setting j\ — —ki + k 2 + k 3 
etc. We choose the current change of variables so that we can use some of the 
computations in [S]. In fact, the set 

(4.14) M* n := fjei: -An < j v - j„ < An, 1 < v, (i < 4} 

is used in [9] . The main result of this subsection is a compact formula for the partial 
sum 

(4.15) D n (x) := e k(x) = e j(t) = : ^n(t), ej(t) := e^*, 

where x and t are related by (14. lip and the middle equality follows from the fact 
that A^* = G„ under this change of variables. In fact, by (|4.11|) and (|4.13p . we 
have 

k-x = kx{t 2 + t 3 ) + k 2 (ti + t 3 ) + k 3 (h + t 2 ) 
= {k 2 + k 3 )h + (h + k 3 )t 2 + (ki + k 2 )t 3 

= j [0'i _ k)h + (h ~ k)h + (j3 - 34,)t 3 ] = -j ■ t 

where in the last step we have used the fact that t 6 R^. The compact formula of 
D n (t) is an essential part of the compact formula for the interpolation function. 

Theorem 4.6. For n > 1, 

D*(t) = e„ +1 (t) - e n (t) - (e° dd (t) - e°^ 2 (t)) , 

where 

e n (t) = f[ s ^ m "' 



sin irti 

2 — 1 



and for n > 1, 



4 . / , „x , 4 



and 



2—1 7—1 V 



erw=fr sin(n+1) ^ y sin S n+3Wj \ ^=0^. 

re w 11 sin 27^ ^ sinfn + lW*,- 

i=l 1 = 1 v ' J 
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Proof. By definition, G n is a subset of H* that contains elements with all indices 
being even integers. For technical reasons, it turns out to be easier to work with 
H* \ G n . In fact, the sum over H* has already been worked out in [S], which is 

e aw = n sin(n+ ; Ktt - n ^ - @ ^ - e «w- 

^— ' - L - L sm7r£,' x± suntU 

j e H* i=l 1 i=l 1 

Thus, we need to find only the sum over odd indices, that is, the sum 
A° dd (t):= £ ej(t), H° dd := H* \ G„. 

jgjflodd 

Just as in [9] , the index set H° dd can be partitioned into four congruent parts, each 
within a parallelepiped, defined by 

K k) ■= {J G Hf : < ji - Jk < An, I G N 4 } 

for k G N4. Furthermore, for each index set J, C J C N4, define 

:= {k G H° dd : k % = kj, Vi, j € J; and < ki - kj < An, V? € J, Vi G N 4 \ J} . 

Then we have 

Using the inclusion-exclusion relation of subsets, we have 

K dd (t)= E (-i) |J|+1 E 

0C,/CN 4 k£H;{ 

Fix j G J, using the fact that t j = — Si^j > we have 
Since k G implies, in particular, ki = kj mod 4, we obtain 

E e * kt = n e e* n E 

keH^ ieN4\J0<fei-fc3<4n ZeN 4 \ J CKfci <n 

keH^ |fc|jodd 

where := X);gn 4 \j^- The ^ as ^ equation needs a few words of explanation: 
if Ak[ = ki — kj, then using the fact that ki = kj, Vi, j G J for k G and 
ki + ki + + &4 = 0, we see that j Y^ieFSAji^l — = which is odd by the 
definition of H^; on the other hand, assume that J2i&i 4 \j H 1S °dd, then we define 
kj = — X)ieN 4 \j K f° r au i G J and define fc; = 4fcZ' + fcj, so that all components of 
k are odd and k G H^. 

The condition that \k\ j is an odd integer means that the last term is not a simple 
product of sums. Setting 

A?W:= E e2 

j—0,j odd 



27T2 Mz 

5 



^■t _ e 2 ™'(l-e 4 "W') 



« 1 _ 4wil2±Z]t 
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we see that, up to a permutation, only products D^D^D^ and D°D E D E are 
possible for triple products (|J| = 3), only D°D E is possible for double products 
(| J | = 2), only D% is possible (|J| = 1), and there is a constant term. Thus, using 
the fact that abc — (a— 1)(6— l)(c — 1) = ab+ ac+ be— a — b — c+ 1, we conclude 
that 

(il,-t2,-t3)6N4 

+ D°(ti) [D%(t 2 )DZ(t 3 )DZ(t 4 ) - (D%(t 2 ) - l)(D*(t 3 ) - 1){D%{U) - 1)] 

+ D°{t 2 ) [D%(ti)D%{ts)D%{U) - (D*(ti) l)(X>*(t 3 ) 1)(D*(U) 1)] 

+ D°(i 8 ) [D^D^D^U) - (D E (h) 1)(^(* 2 ) - 1)(^(*4) - 1)] 

+ D£(f 4 ) [£ r f (t!)Df (* 2 )^(t 3 ) - (^(ti) - l)(Af (h) l)(D E (t 3 ) ~ 1)] , 

where the first sum is over all distinct triple integers in N4. 

Assume that n is an even integer. A quick computation shows that 

n°(t \n E (t \r> E (+ \n E (t \ - 17 sin7T ( n + 2 ) t * simrnh 
D n ( tl )D n (t 2 )D n (t 2 )D n (t 4 ) - [I s 

J=2 

Furthermore, we see that 

D%(t 2 )D%(t 3 )D%(U) - D%(h)(D E (t 2 ) l)(D E (h) 1)(D E (U) 1) 
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n siriTTTifj i7rnti _ e - etTi-mi j-r 

sin27rf, i_ e 4«*i 11 



i=2 



eiirnti 



4 



sin7rn£i sin7r(n — 2)ti 



sin 27rf 7 - sin 27rfi 

j=2 



Adding the two terms together and then summing over the permutation of the 
sum, we end up the formula for D° dd (t) when n is even. The case of n odd can be 
handled similarly. □ 

Let us write down explicitly the function D n (x) defined in (14.15|) in x- variables. 
Using the elementary trigonometric identity and (|4.11|) . we see that 

4 

4 JJ^ saxemti = (cos ock(x 2 — X\) — cos a7TX3)(cos an{x 2 + x±) — cos 0:71x3) 



i=l 



cos 2 axi + cos 2 ax 2 + cos 2 ax 3 — 2 cos axi cos ax 2 cos ax 3 — 1, 



so that we end up with the compact formula 

(4.16) D n (x) = & n+1 (x) - e B (x) - (e° dd (x) - e°i d 2 (x)) , 



where 

e n (x) 

and 



cos nirxi + cos nirx 2 + cos nirx 3 — 2 cos wnxi cos n7TX2 cos n7rx3 — 1 
cos 2 7rai + cos 2 irx 2 + cos 2 77x3 — 2 cos 7txi cos ttx 2 cos 7TX3 — 1 

9° dd (x) = 6„ +2 (2x) V . f 1 "^ - , if n = even, 

.7=1 J 



9° dd (t ) ^9. +1 (2x)V Sm {" + ^ , ifn = odd, 

.7=1 J 
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in which U is given in terms of Xj in (|4.11|) . As a result of this explicit expression, 
we see that D n {x) is an even function in each x\. 

4.4. Boundary of the rhombic dodecahedron. In order to develop the inter- 
polation on the set X* , we will need to understand the structure of the points on 
the boundary of = Z 3 D £Ib- As J7b is a rhombic dodecahedron, we need to 
understand the boundary of this 12-face polyhedron, which has been studied in 
detail in [9] . In this subsection, we state the necessary definitions and notations on 
the boundary of fls, so that the exposition is self-contained. We refer to further 
details and proofs to [S]. 

Again we use homogeneous coordinates. For i,j G N4 := {1,2,3,4} and i 7^ j, 
the (closed) faces of £Ib are 

F id ={te~n H : U-tj = 1}. 

There are a total 2( 2 ) = 12 distinct Fij, each represents one face of the rhombic 
dodecahedron. For nonempty subsets J, J of N4, define 

^/,J := P| F id = {t G n H ■ tj =ti-l, for all i £ I,j e J} . 

It is shown in [9] that Sl^j = if and only if ID J ^ 0, and n ^/ 2 ,j 2 = if 

I1UI2 — I and J1UJ2 = J. These sets describe the intersections of faces, which can 
then be used to describe the edges, which are intersections of faces, and vertices, 
which are intersections of edges. Let 

JC := {(I, J) : I, J C N 4 ; I n J = 0} , 

/C := {(I, J) G JC : i < 3, for all £ {I, J) } . 

We now define, for each (I, J) G JC, the boundary element B/.,/ of the dodecahedron, 

B 7iJ := {t G : t g ^ 7l .. 7l for all (J x , J x ) G JC with |/| + \J\ < \h\ + | Ji|} ; 

it is called a face if |/| + | J\ = 2, an edge if |7| + \ J\ = 3, and a vertex if |/| + \ J\ = 4. 
By definition, the elements for faces and edges are without boundary, which implies 
that Bj^j n Mi^ji = if I ^ I\ and J ^ J\. In particular, it follows that B{i},{j} = 
i 7 !^ and, for example, ®{i},{j,fc} — {Fi,j H -Fi.fc) for distinct integers i, j, fc G N4. 

Let G = S4, denote the permutation group of four elements and let denote 
the element in Q that interchanges i and j; then tcry = t — (ti — tj)e^j. For a 
nonempty set / C N4, define Gi := {aij : i,j £ J}, where we take <jjj = aji and 
take (jjj as the identity element. It follows that Gi forms a subgroup of G of order 
|/|. For (I, J) G /C, we then define 

(4.17) [Bj,j] := |J Bi,j(r. 

It turns out that [Bj,j] consists of exactly those boundary elements that can be 
obtained from B/ ; j by congruent modulus B, and [V>i,.j] n [B^.jJ — if (I, J) ^ 
{I\,Ji) for (I, J) G /Co and (/1, Ji) G /Co- More importantly, we define, for < 
i,j < i + 3 < 4, 

(4 18) ® ij := IJ t B '^ with : = ^) e /Co : |/| = i, |J| = j} • 

(/,./)e/c^ 
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Then the boundary of Qb can be decomposed as 

Tl H \n° H = |J B 7 , 7 = |J B*'- 7 ' . 

(7,J)£/C 0<i,j<i+j<4 

The main complication is the case of |/| + |J| = 2, for which we have, for example, 

(4-19) PB{1},{2,3}] = B {1},{2,3} U B {2},{1,3} U B {3},{1,2}- 

The other cases can be written down similarly. Furthermore, we have 

B {1} : {2,4} = ®{1},{2 : 3} <J 34, ®{1,2},{4} = ^{l^},^} ^) 

(4.20) B{1},{3,4} = B{i},{ 2 ,3}°"24, 5{1,3},{4} = B{1,2},{3}0'230'34, 

B {2},{3,4} = ®{l,2},{3} cr 120'24, B{2,3},{4} = B{1,2},{3}°'130'34, 

with 

2 B{i},{2,3} = {(*, « - 1, t - 1, 2 - 3t) : | < t < f } , 

B { i, 2}i{3} = {(l-t,l-t,-t,3t-2): £<*<!}• 
If |J| + | J | = 2 then B 7; j = B^j^jj is a face and 

B 1 ' 1 = [B{i},{2 } ] U [B { i } , {3} ] U [B { i}, {4} ] U [B { 2}, { 3 } ] U [B {2} , {4} ] U [B {3} , {4} ] 
If |/| + | J| = 3 then B/ ; j is an edge and we have 

^ B 1 ' 2 = [B {1}i{2 ,3}] U [B {1} , {2i 4}] U [B{i},{ 3 , 4 }] U pB{2},{3,4}], 

B 2,1 = [B{ 1: 2},{3}] U [B{1,2},{4}] U [B{ li3 } : {4}] U [B{2,3},{4}]- 

If |/| + |J| =4, then 

Recall that G„ is A£* = Z 3 n Ob in homogeneous coordinates. We now consider 
the decomposition of the boundary of G„ according to the boundary elements of 
the rhombic dodecahedron. First we denote by G° the points inside G n , 

G° := {j e G„ : -An < j v - j„ < An, 1 < v, n < 4} = {j e G„ : ± e O s } . 

We further define, for < i, j < i + j < 4, 

(4.24) G«:={kgG B :ieB«} 

The set G^- 7 describes those points j in G„ such that ^ are in B 1 ^ of c3Sl s . It is 
easy to see that G£ j n G% 1 = if i ^ k, j ^ I and 

|J G^' = G„\G°. 

0<i,j<i+j<4 

4.5. Interpolation by trigonometric polynomials. We first apply the general 
theory from Section 2 to our set up with as a rhombic dodecahedron. 

Theorem 4.7. For n > 1 define 

(4.25) I n f(x) := ~ £)- M*) : = ^ E ^C*)- 
TTien /or eac/i j e X„, I n {j^) = /(^)- 
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Proof. By @U, = /(£) for j g X n is equivalent to I n f(B~ tr l) = f(B~ tr l) 

for Z g A„. Moreover, /„/ can be rewritten as 

Mix) = E f{B- tr m n {x-B- tr j). 

ieA„ 

Hence, this theorem is a special case of Theorem 12.21 □ 



Next we consider interpolation on the symmetric set of points X*. For this we 
need to modify the kernel function Recall that, under the change of variables 
(I4.13|) , A],* becomes G„ in homogeneous coordinates. We define 



4> 



X*) ■= i E /$°*(*) = ± E ^ B) ej(t), 



2n 3 

^eAl* J - 



where x and t are related by (|4.1ip . /i^ is defined by li^ under the change of 



indices (|335|) . and zxj nJ = 1 if j g G°, ixj nJ = if j g G^'; more explicitly 



,(») 



1, 

j € GM, 
jgG^UG^, 
jgG^UG 3 /, 



^ 6 

For each k on the boundary of X*, that is, ^- on the boundary of [— i, i] 3 , let 

(4.26) 5 k :={jei; : iEi modZ 3 }, 

which contains the points on the boundary of X* that are congruent to k under 
integer translations. 

Theorem 4.8. For n > 1 define 

(4.27) ^/(x) := E f(&)Rk(x), R k (x) := - A). 
TTien /or eac/i j g X*, 

(4.28) J*/(i) 



/ ( 2n ) ' 



fees, 

In homogeneous coordinates, the function <&* n {x) — 3>£(t) is a real function and it 
satisfies 



Kit) = 



(4.29) 



1 

in 3 



\ [Kit) + D^it)) - \ E ^Sr^ E ^2,{nt + Lf jt.) 



3 ^— ' sin 27rti/ 

!/=l j = l 



1 „ 1 cos27mi,-, if 

- E cos2wnit fl + U)--{^- 1 31 J 

3 i<^t<4 2 1° l f 



n even 
n odd 
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from which the formula for <I>* (x) follows from ([4. lip and (|4. 16|> . 

Proof. By (|4.1[) . we need to verify the interpolation at the points B~ tr l for I E A*. 
By definition, we can write 

fifc(B- tr 4l ^ Il) e,(i?- tr (/ - fc)). 

It is easy to see that v tr B~ tr l = jt(jih + hh + jzh) if v is related to j by (|4. 13|) . 
Hence, as in the proof of Theorem 3.15 in [9], we conclude that 

Rk(B- tr l) = ^ e,(B- tr a-fc)), 

Now, for I, k £ A*, there exist p S A„ and g S Z 3 such that I — k = p ± B tr q. 
Consequently, it follows from (|2.6p that 



Rk(B- tr l) = ^3 E ^C 5 "^) = <W 



By (|4.ip , we have verified that 
(4-30) R k {±) = 



1, i- = #- mod Z 3 

' 2n 2n 

0, otherwise, 



which proves the interpolation part of the theorem. 

In order to prove the compact formula, we start with the following formula that 
can be established exactly as in the proof of Theorem 3.15 in [9]: 



\{K$) + K-xW>)-\ E ^(t) 
-\ E Mt)-\ E M*) 

fceG^UG 3 / fcGG 2 ' 2 

Let us define G£ J := {k e G n : ^ e B 7 ,j} for /,JcN 4 and also define [G^ J ] := 
{k € G„ : i € [Bj,j]}. It follows from (j4T8)) . and (|4~24|) that 

G«= |J [G'/] and [G£ J ] = (J G£ J <r. 

In order to compute the sums in (|4.3ip , we need to use the detail description of the 
boundary elements of fig in the previous subsection. The computation is parallel 
to the proof of Theorem 3.15 in [S], in which the similar computation with G„ 
replaced by H„ is carried out. Thus, we shall be brief. 

Using t\ + ti + £3 + ti = and the explicit description of B^ 1 ^ 2 ' 3 }, we get 



v Mt)= E E e 

keGi 1} ' {2 ' 3} keGi 2} - {1 - 3} k G Gi 3 >-- 



e -2irijU f e 2nni(ti+U) _|_ e 2n7ri(t 2 +*4) _)_ e 2n7ri(i 3 +t 4 )^ 

1 j even 

l27r L^~J^4 c ^ 2 TT^| S+l|t 4 / e 27rm(t 1 +t 4 ) + e 27rm(i 2 +t 4 ) + e 27rm(t 3 +t 4 )\ 

sin 2-7rf4 V / ' 
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Similarly, we also have 

E Mt)= E e¥kt + E e¥kt + E e¥kt 

kelG* 1 ' 2 ^*^] keGi 1} ' {2 - 3} keGi 2 >> {1 ' 3 > keG< 3 >' {I > 2} 



sin27r L^jt4 e27r . 



L^JU ^ e -27rm(ti+t 4 ) _|_ e -27rjn(t2+t 4 ) _|_ e -27rjn(t 3 +t 4 )~j 



sin 27r<4 

From these and their permutations, we can compute the sum over G*' 2 and G 2,1 . 
Putting them together, we obtain 

E Mt) = 2± ^g^Ecos2.K. + Lf 

keG^UG 2 ' 1 v=x " i = 1 

Using (|4.23p , we see that, G 2,2 = {(2n, 2n, —2n, —2n)a : <r G 5} and, if n is even 
then G^' 3 = {(n,n,n, —3n)a : cr E Q} and G 3,1 = {(3n, — n, — n, — n)a : a e C*}, 
whereas if n is odd, then G^ 3 = G 3 ' 1 = 0. As a result, it follows that 

E tfk(t)= E e 2 ™^ +t ") = E cos27rn(V + ^), 

keG 2,2 l<fi<i/<4 l</i<i/<4 

where we have used the fact that t\ + t 2 + £3 + £4 = 0, and 

4 4 
E <t>k (t) = E (e 27rmtj + cr 2 ™ 1 *; ) = 2 E cos 27rn*j , 

keG^UG 3-1 3 = 1 3 = 1 

if n is even, whereas it is equal to if n is odd. 

Putting all these into (|4.31[) completes the proof. □ 

Theorem 4.9. Let ||/*||oc denote the norm of the operator I* : C([— 5, ^] 3 ) 
C([— |, ^] 3 )- T/ien ttere is a constant c, independent of n, such that 

3 



II ^ II 00 < c(logn) 
Proof. Following the standard procedure, we see that 



\I*\\oo= max E l $ n(*- &)|. 



3:61 2 ' 2J feeX » 



Using the formula of <£>* in (14. 29[) . it is easy to see that it suffices to prove that 

max \D*(x - < c(logn) 3 , n > 0. 

tl 2,21 feex* 

Furthermore, using the explicit formula of D° dd (t) and (3.19) in 0, we see that 
our main task is to estimate the sums in the form of 



Ifi 2 3} : = tt^7 max 



2n 3 teg 



sinTrnfr - *$a) sin7rn(t 2 - ^) sin7rn(t 3 - 



sm7r(ti - ^)sin^(t 2 - ^)sin^ 3 - ^ 



2n 



and three other similar estimates Zn 2,4}j ^{1,3,4} and i/2,3,4}, respectively, as well 
as similar sums in which the denominator becomes product of sin27r(ii — j 1 ) and 
n in the numerator is replace by n + 1 or n + 2. Here Q is the image of [—1, l] 3 
under the mapping (I4.11[k that is, 

Q = {t G M% : -\ < t 1 +t 2 ,t 2 +t 3 ,t 3 +t 1 < |}. 



25 



Changing the summation indices and enlarging the set X*, we see that 

\ 3 

sinn7r(t - ^) 

t6[-l,l] ' 



/ 1 2n 

Ifi 2 3\ < 4 max — > 



sin7r ( t _ *) 



< c(logn) 3 , 



where the last step follows from the standard estimate of one variable (cf. [301 Vol. 
II, p. 19]). □ 

4.6. Interpolation by algebraic polynomials. The main outcome of Theorem 
14.71 in the previous section is that we can derive a genuine interpolation by trigono- 
metric polynomials based on the set of points in : k € S n } defined at (|4.6[) . 
The development below is similar to the case of d — 2. We define 

£6{-l,l} 3 

Theorem 4.10. For n > define 

fc6S„ 

iw'i/l given in (|4. 9[) . TTiera £ n / E T n is even in each of its variables and it 
satisfies 

Proof. As shown in (|4.30|) . Rk{x) := ® n (x- satisfies Rk{4~) = 1 when fc = j 
mod 2?iZ 3 and otherwise. Hence, if j € S° then {VRk){j^) = ^Rk(^) = 

[^i™^ _1 ^fe,i- If J <= \ 3°, then we need to consider several cases, depending on 
how many components of j are zero, which determines how many distinct terms 
are in the sum (VRk)(-S~) and how many distinct k can be obtained from j by 
congruent in Z 3 . For example, if j E and none of the components of j are zero, 
then there are 2 elements in Sj, j and the one in the opposite face, and the sum 

FRk(^) contains 8 terms, so that (PRk)(J-) = \5j_k — ] -1 ^fe,j- The other 
cases can be verified similarly, just as in the case of d = 2. We omit the details. □ 

The above theorem yields immediately interpolation by algebraic polynomials 
upon applying the change of variables (|4.5|) . Recall r„ defined in l|4.7j) and the 
polynomial subspace 

II* = span{s^S2 2 S3 3 : h,k 2 ,k 3 > 0, h + kj <n,l< i,j < 3}. 
Theorem 4.11. For n>0, let 

£nf{s)= Y f(z k )t k (s), e%(s)=ik(x) with s = cos2?rx. 

z k er n 

Then C n f E II* and it satisfies C n f{zk) — f{zk) for all £ r„. 

This theorem follows immediately from the change of variables (|4.5[) . The explicit 
compact formula of £k(x), thus £%(s), can be derived from Theorem 14.81 

The theorem states that the interpolation space for the point set r„ is exactly 
II* , which consists of monomials that have indices in the positive quadrant of the 
rhombic dodecahedron, as depicted in Figure 3 below. 

The set T„ consists of roughly n 3 /4(l + C(n _1 ) points. The interpolation poly- 
nomial C n f E n* is about a total degree of 3n/2. The compact formula of the 
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(0,0, n) 




Figure 3. Index set of II* 



fundamental interpolation polynomial provides a convenient way of evaluating the 
interpolation polynomial. Furthermore, the Lebesgue constant of this interpolation 
process remains at the order of (logn) 3 , as the consequence of Theorem 14.91 and 
the change of variables. 

Corollary 4.12. Let |j£„||oo denote the operator norm of C n : C([— 1,1] ) >— > 
C([— 1,1] ). Then there is a constant c, independent of n, such that 

||Ai||oo < c(logn) 3 . 
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